Four-well tunneling states and elastic response of clathrates 
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We present resonant ultrasound elastic constant measurements of the clathrate compounds 
EugGai6Ge3o and SrsGaieGeao. The elastic response of the Eu clathrate provides clear evidence 
for the existence of a new type of four-well tunneling states, described by two nearly degenerate 
four level systems (FLS). The FLS's are closely linked with the fourfold split positions of Eu known 
from neutron diffraction density profiles. Using a realistic potential we estimate the tunneling fre- 
quencies and show that the energy gap between the two FLS's is of the same order as the Einstein 
oscillator frequency. This explains why the observed harmonic oscillator type specific heat is not 
modified by tunneling states. In addition the quadrupolar interaction of FLS's with elastic strains 
explains the pronounced depression observed in elastic constant measurements. In the case of the 
Sr clathrate, we show that the shallow dip in the elastic constant C44 is explained using the same 
type of quadrupolar interaction with a soft Einstein mode instead of a FLS. 

PACS numbers: 62.30.+d,65.40.-b,66.35.+a 
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Clathrate compounds are composed of polyhedral 
cages with 12-16 faces formed by Si, Ge or Ga atoms 0. 
The oversized cages accommodate guest atoms like Ba, 
Sr, or Eu which exhibit large amplitude anharmonic 
motion. This leads to a strong scattering of acoustic 
phonons and a small thermal conductivity. In many cases 
the electrical conductivity remains high and the result- 
ing large thermoelectric figure of merit has led to sus- 
tained interest in this class of compounds. The glass- like 
low temperature thermal conductivity 0, and ultra- 
sound attenuation 0, |j| have been interpreted as evi- 
dence for the existence of tunneling states. These are 
commonly discussed within a generic two well poten- 
tial model leading to two level system (TLS). The ther- 
mal conductivity is more glass-like for clathrates with 
larger atomic displacement parameters (ADP) of the 
guest atoms. For the clathrates discussed here the ADP's 
are very large. The nuclear density map for Eu(2) in 
EusGaigGeao clearly shows the distinct symmetrically 
related four maxima, away from the cage center 0, [3, 
leading to the fourfold split site model. Large static 
ADP for Sr(2), and considerably smaller for Ba(2), in 
Sr§Gai6Ge3o and BasGaifiGeao, respectively, point also 
to the split site model 0, |g, H| , although the nuclear den- 
sity maps do not show distinct fourfold maxima. 

In this letter we present elastic constant measure- 
ments on XgGai6Ge3o (X = Eu, Sr) that show anomalies 
that cannot be naturally explained within a generic TLS 
model. Instead we propose two new models that we feel 
are physically appealing because they draw their inspira- 
tion from the nuclear density map results. 

For EusGaigGeao we introduce a four- well potential 
model with the proper symmetry of split sites that de- 
scribes a rotor type tunneling motion. We show that for 
appropriate potential parameters two partially degener- 



ate FLS's exist. This model explains both the specific 
heat and elastic anomalies and is compatible with neu- 
tron diffraction. 

For SrgGai6Ge3o the elastic response is explained 
using the same quadrupolar interaction with elastic 
strains, but with local Sr vibrations. This model may 
be also considered in other compounds with soft local 
modes, such as skutterudites, where a similar shallow 
dip in the elastic constant is observed 0, without 
pronounced feature in the specific heat. 

The elastic constants reported in this work were mea- 
sured using Resonant Ultrasound Spectroscopy (RUS) 
[Io| . RUS is a novel technique for determining the elastic 
moduli of solids, based on the measurement of the reso- 
nances of a freely vibrating body: in a RUS experiment, 
the mechanical resonances of a freely vibrating solid of 
known shape are measured, and an iteration procedure 
is used to 'match' the measured lines with the calculated 
spectrum. This allows determination of all elastic con- 
stants of the solid from a single frequency scan, eliminat- 
ing the need for separate measurements to probe different 
moduli. RUS measurements were taken as a function of 
temperature using a homemade probe that fits in a Phys- 
ical Property Measurement System (PPMS) from Quan- 
tum Design. Heat capacity measurements were also per- 
formed in the PPMS using the apparatus from Quantum 
Design. 

The four well potential is modeled, in accord with nu- 
clear density data, as 
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The harmonic oscillator (HO) potential of the last term 
is modified by the four-well potential term. Vo is the 
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FIG. 1: Left: evolution of p with potential barrier Vo- Dashed 
lines represent nondegenerate (for Vo > 0) and the solid lines 
doubly degenerate levels. The formation of lower FLSi and 
upper FLS2, may be followed. Right: eigenvalues \ nm and 
potential V(0,p) on the circle p — po with parameters ap- 
propriate for EusGaieGeao. Lower FLSi: fat line, doublets: 
thick line, nondegenerate state: thin line. Energies in units of 
K. For every eigenvalue p m , there is a sequence of equidistant 
eigenvalues \ n ,m of the radial equation (Eq0 separated by 
4 x 9/2 = 60 K. 
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Sr [11] 


87.6 


0.35 


50 


0.11 


2.23 


56 


Ba [3] 


137.3 


0.20 


70 


2.19 


4.45 


110 


Eu [3,6, nj 


152.0 


0.43 


30 


0.05 


0.85 


22 



TABLE I: The parameters of "guest" atoms for the three 
clathrates. A is the atomic mass number. Ai, A2, and A3 
are the tunneling energies of the FLSi, FLS2, and the gap 
between the two, respectively, obtained for Vo = 30. pp for 
Eu is the average from the three references. (In Ref. .3] the 
position for Ba(2) is incorrect. Here the correct po is used.) 



coordinate and A 
eigenvalues 
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2E/(ftB<d) denotes dimensionless 
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The maxima of the radial wave functions correspond 
to the nuclear density maxima. The tunneling splitting 
Ai in units of Kelvin for the lowest FLSi is given by 



Ai 



potential barrier in the azimuthal direction. The Vi 
term determines the positions of the minima, located in 
the directions <f> — (2n + 1)tt/4, for n = {0,1,2,3}, at 
Pmin — \j2V\jK. It also represents the infinite poten- 
tial barrier across the center of the cage. This means 
that there is no tunneling through (and no nuclear den- 
sity in) the center, consistent with the nuclear density 
map for EusGai6Ge3o [EQ- The HO force constant is 
given by Einstein temperature, 0, of the soft local mode 
of the guest atom of mass M, K = M (ksQ/fi) 2 ■ may 
be determined from ADP, Raman, or specific heat data. 
The wave function may be written, after separation of 
the variables, as the product of angular and radial parts, 
ip(p, 4>)n,m = R(p)n, m ^((f>)m- The former obeys 
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where Vo 

(0) 

ues fXm = 



Tpr^O is dimensionless. For Vo = eigenval- 
m 2 , eigenfunctions $^(0) = cos mip, sinmcj) 
and angular momentum mh correspond to the free rotor. 
This has a formal analogy to directional rotors formed 
by guest molecules in clathrates 01 • Introduction of Vo 
breaks the full azimuthal symmetry to fourfold symme- 
try, lifts the twofold degeneracy of states with even values 
of m and leads to the formation of FLS with increasing 
Vo, shown in Fig.^ The radial wave function satisfies 
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Here, x 2 = (M K /h 2 ) 1 / 2 p 2 is the dimensionless radial 



Here, po is the split site radial distance from the center, 
corresponding to the maximum of the n — radial wave 
function. The tunneling frequency A2 of the upper FLS2 
and its energy relative to the lower one, A3, may be 
calculated using Q}. A3 is a measure for the size of 
the azimuthal corrugation in the radial minimum. The 
low lying energy levels for the case of Eu clathrate, are 
shown in Fig.^ The parameters for the three clathrates 
are summarized in Table [I] We note that the nuclear 
density profiles for Sr and Ba do not correspond to the 
FLS wave functions presented in the right part of Fig. 
Our new model of azimuthal four-well tunneling enforces 
a re-interpretation of low temperature properties in Eu 
clathrate. 

Using the eigenvalues from ©, the specific heat, C(T), 
may be numerically calculated. For the free rotor (Vo = 
0) closely spaced levels (Fig. lead to a constant C(T) 
for moderately low T (Eq. [3] and Fig. • Turning on the 
four well term, Vo, leads to the "bunching" of the levels 
into two FLS's, with large energy gap A3 between them. 
Since A3 ~ O for Eu potential parameters and Ai is very 
small, C(T) for the pure HO and FLS model deviate very 
little. They are equally acceptable, in contrast to the 
conventional TLS model. In addition, the contribution 
of guest atoms to the total C(T) at high temperatures 
is only around 10%, so that the deviations in the high 
temperature region are negligible, as seen in Fig. [2] 

Now we show that our new FLS model in addition leads 
to characteristic elastic constant depressions observed in 
EusGaieGeso by RUS measurements which the pure HO 
model cannot explain. Assuming a quadrupolar inter- 
action V cos(2<p) of FLS's with the lattice, the matrix 
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FIG. 2: Top left: specific heat contribution C v of a guest 
atom. Dashed line: unperturbed 2D-HO, dash-dotted line: 
2D free rotor with Vi = 307 and Vo = 0, solid line: FLS's 
with Vi — 292 and Vo — 30. These parameters are in accord 
with the case of Eu (TableQJ. Top right: contour plots of the 
ground state wavefunctions for FLSi and FLS2. The shaded 
background corresponds to the four well potential. Bottom: 
measured specific heat (circles) of EusGai6Ge3o from Ref. 
Solid line is the calculated total specific heat Cv including the 
sum of 6 FLS/cell, HO contributions and the Debye term. It 
is almost identical to the fit using the pure Einstein oscillator 
for the guest atom Eu. The peak around 30 K is due to 
the ferromagnetic phase transition. Inset: schematic repre- 
sentation of FLSi,2- Arrows indicate transitions within each 
FLS system due to the quadrupolar interaction. 



elements between the eigenstates of Eq. can be nu- 
merically calculated. The interaction matrix is nearly of 
block diagonal form with negligible off-diagonal terms, 
which is represented as direct sum of the two FLS's: 



( liX 

^ l2 x 
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with X — a x ® a x (a x — Pauli matrix) and 71, 72 the 
quadrupolar coupling strengths for FLSi and FLS2, re- 
spectively (71 > 72). The possible transitions are shown 
in the inset of Fig. [21 The interaction term in Eq. (JSJ) 
is coupled to the elastic strain e p leading to the total 
Hamiltonian 



H = lvc° u el 
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pie volume and Hfls the noninteracting two-FLS Hamil- 
tonian. The renormalised elastic constant is then given 
by the static quadrupolar susceptibility of the FLS's ac- 
cording to Ac M = (N c /V )xq where N c — 6 is the num- 
ber of Eu(2) atoms per unit cell of volume Vq and 



XQ 



E 



(m\V mt \n}\^ 



(Pm ~ Pn) 



^0,n ~ Arj,r, 

^{ E \(™\V m t\n)\ 2 p m -(V mt ) 2 }(l) 



k B T 



Here p m is the occupation factor for the FLSi, 2 tunnel- 
ing state I to). The first (Van Vleck) contribution is due 
to inelastic transitions within FLSi and FLS2. The sec- 
ond (Curie) term is due to elastic transitions within the 
degenerate doublets of both FLSi and FLS2. 

The FLS's interact with the phonon bath. For a suffi- 
ciently strong interaction compared to the tunneling en- 
ergy the incoherent regime is entered. The crossover tem- 
perature is given by T* = ^\J^\Zpv^^ Ai is 
an order of magnitude smaller than A2 (Table and 
correspondingly, T* for the FLSi is significantly lower. 
In the coherent regime T< T* the change of \Q is d ue 
to the relative change of potential minima between the 
four wells. Because A x is small its coherent tunneling 
states are destroyed by the coupling to the phonon bath 
already at T ~ IK. Then the guest atom may be consid- 
ered as statistically distributed over the four wells reduc- 
ing the elastic response to that of four independent wells. 
This is zero if only the constant shift of a single well is 
considered. In contrast, for Eu parameters FLS2 is on 
the boundary of the coherent regime in the temperature 
range of the depression. Therefore only FLS2 contributes 
to the quadrupolar response, explicitly: 
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Here, c° is the background elastic constant, V is the sam- 



The depression in c M is caused by the thermal depopula- 
tion of the relaxational FLS2 , described by the exponen- 
tial factor in front of Eq. (JHJ. Within a conventional TLS 
model it may only be described by assuming an unrea- 
sonably large tunneling frequency, incompatible with the 
specific heat data, and unequal equlibrium well depths 
which in the present context is unphysical because it vi- 
olets the fourfold symmetry. 

The results for the Eu clathrate are shown in Fig. 
Our fit determines 7^2 ~ 0.073 eV, with A 3 /k B = 22 K. 
This value is used above to estimate T* and to jus- 
tify the neglect of the FLSi contribution. The depres- 
sion of c M is a direct signature of the quadrupolar cou- 
pling to relaxational processes within the upper FLS2. 
Contrary to TLS model our FLS tunneling model ex- 
plains both specific heat and elastic constant anomalies 
in EugGaigGeao without contradiction and is consistent 
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FIG. 3: The fit of the elastic constants of Eu clathrate in- 
cluding interaction with FLS2 (solid thick line). Background 
variation corresponds to Varshni function [T^| (dotted line). 
The depression with respect to the Varshni background is due 
to quadrupolar coupling to the upper FLS2. Note the differ- 
ent scales for c M . 

with the nuclear density profile from the neutron diffrac- 
tion. This presents strong evidence that a new type of 
FLS tunneling states due to azimuthal tunneling of guest 
atoms has been discovered in Eu clathrate. 

In the case of Sr 8 Gai6Ge3o the split sites are less dis- 
tinct than in EusGaifiGeso, with much of the nuclear 
density in the center |6| , and the four- well potential pre- 
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FIG. 4: Elastic constant C44 in SrgGai6Ge3o. Fitted curve 
(solid thick line) is obtained using the Varshni function for the 
background (dotted line) and the model of the quadrupole in- 
teraction of the local Einstein mode with the elastic strain ex- 
plained in the text. The interaction constant jho = 0.015 eV. 
Inset: schematic representation of the transitions between the 
states of 2D-HO due to the quadrupolar interaction. Upper 
inset: the contour plot of the quadrupolar interaction poten- 
tial Vf n ?=lHOP 2 sin(2</>) exp [~p 2 /(2og)], a = hj^KM. 



sented here is not appropriate. Due to the large overlap, 
as may be deduced from the nuclear density profile, the 
tunneling frequency may be very large. The split-site 
model actually has the same agreement with the data 
as the single site model 0. Therefore, we propose a 
different model for Sr 8 Gai 6 Ge3o. We assume that the Sr 
atoms sit in the shallow HO potential that produces a soft 
Einstein mode. However we now propose a quadrupolar 
form of interaction which rapidly decays at larger dis- 
tances from the center. This is a crucial condition for a 
small dip in the elastic constant to appear due to tran- 
sitions between the excited degenerate states of 2D-HO. 
In the upper inset of Fig. 0] this interaction potential is 
shown. Although the four well guest atom potential for 
Sr has not yet developed as in the Eu tendency 
is present in the model quadrupolar interaction potential 
which leads to a satisfactory agreement with RUS data 
in Fig.H 
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